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Abstract 



We generalize our previous model [9] to an 0(N) symmetric two- 
dimensional model which possesses chiral symmetry breaking {{ipip) conden- 
sate) and superconducting (Cooper pair {ijj-tjj) condensates) phases at large-N. 
At zero temperature and density, the model can be solved analytically in the 
large-N limit. We perform the renormalization explicitly and obtain a closed 
form expression of the effective potential. There exists a renormalization 
group invariant parameter S that determines which of the {'ip'ip) {S > 0) or 
{tp'<p) < 0) condensates exist in the vacuum. At finite temperatures and 
densities, we map out the phase structure of the model by a detailed numerical 
analysis of the renormalized effective potential. For 6 positive and sufficiently 
large, the phase diagram in the /x-T (chemical potential-temperature) plane 
exactly mimics the features expected for QCD with two light flavors of quarks. 
At low temperatures there exists low-/x chiral symmetry breaking and high-/x 
Cooper pair condensate regions which are separated by a first-order phase 
transition. At high /x, when the temperature is raised, the system undergoes 
a second-order phase transition from the superconducting phase to an unbro- 
ken phase in which both condensates vanish. For a range of values of S the 
theory possesses a tricritical point (ptc and Ttc); for fi > fxtc (/U < fitc) the 
phase transition from the low temperature chiral symmetry breaking phase 
to unbroken phase is first-order (second-order). For the range of 6 in which 
the system mimics QCD, we expect the model to be useful for the investiga- 
tion of dynamical aspects of nonequilibrium phase transitions, and to provide 
information relevant to the study of relativistic heavy ion collisions and the 
dense interiors of neutron stars. 



PACS numbers: 
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I. INTRODUCTION 



The phase structure of QCD at non-zero temperature and baryon density is important 
for the physics of neutron stars and heavy ion colhsions. 

An approximate phase structure for QCD with two massless quarks has been mapped out 
in various mean field and perturbative approximations and a rich structure has emerged. For 
a recent review, see . In addition to the well known chiral symmetry broken and restored 
phases recent investigations have revealed the possibility of a color superconducting phase at 
low temperatures and relatively high densities [|]. In the chiral condensation regime 
at zero chemical potential, the phase transition to the unbroken mode is second order as we 
raise the temperature As we increase the chemical potential at fixed low temperature 
there is a first order transition to a superconducting phase. There is also a regime where 
as we increase the temperature the phase transition from the chirally broken phase to the 
unbroken phase is first order, so that somewhere along the line separating these phases there 
is a tricritical point [||. These results are summarized in 

One of the more interesting questions is what happens in a dynamical situation such 
as a heavy ion collision, in which the system traverses the various phase transitions as it 
expands and cools. One would like to study the correlation functions to see whether there is 
qualitatively different behavior in crossing the first order or second order transition regions 
as a function of the proper evolution time and whether this difference would lead to some 
interesting experimental signatures at RHIC. 

In order to get a better handle on this latter question, we propose a simple 1+1 dimen- 
sional model which contains several of the features of two flavor massless QCD in mean field 
approximation (similar phase structure and asymptotic freedom) which would lend itself 
to dynamical computational simulations pertinent to the (one dimensional) expansion of a 
Lorentz contracted disc of quark matter. Thus we hope to explore the behavior of a system 
of "quarks" evolving through first and second order phase transitions to either a final state 
with chiral condensates or to a superconducting final state. This calculation would be sim- 
ilar in spirit to those that were done in exploring the chiral phase transition in the linear 
sigma model 

In this paper we restrict ourselves to mapping out the phase structure of our model at 
finite temperature and chemical potential in a large- approximation, for use in obtaining 
initial conditions for our future dynamical calculations. Our simple model combines the 
Gross- Neveu model with a model for Cooper pairs that we introduced recently ||^. It 
turns out to have many of the features of QCD that we ultimately want to capture in more 
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realistic calculations. Namely, the theory in the Gross-Neveu sector has the same second- 
order, first-order, tricritical point behavior for the chiral condensate < gg > as a function of 
temperature and chemical potential [rU| |n| as QCD with two massless flavors. Adding 
the second interaction also adds a new phase where there is superconductivity at some finite 
chemical potential as also expected in QCD with two massless flavors. The model has a well 
defined expansion and is asymptotically free so that it does not suffer from the cutoff 
dependences of 3 + 1 dimensional effective field theories considered by others ||^ . 

Thus we will be investigating a 1+1 dimensional model governed by two independent 
couplings: the original Gross-Neveu term |^ that promotes the condensation of (gg), and the 
term considered in our earlier paper 0that produces a {qq) condensate. We first determine 
the unrenormalized effective potential at leading order in a large-N expansion, by introducing 
collective coordinates for the gg, and gg operators, and integrating out the fermions in the 



usual fashion by using the Hubbard-Stratonovich trick [12 



At zero temperature and chemical potential we derive a closed-form analytic expression 
for the renormalized effective potential. We find that there is one dimensionless parameter 5, 
independent of the renormalization scale, whose value determines which of the condensates 
is present. This situation might be described as "partial dimensional transmutation": the 
unrenormalized theory has two bare couplings whereas the renormalized one has a renor- 
malization scale, which is arbitrary (but which can be related to the physical fermion or 
Cooper pair gap mass) , and a dimensionless parameter 5 independent of this scale, that 
controls the physics. We find that the gap equations have three types of solution: two in 
which one or the other of the condensates vanish, and a third, mixed case, in which both 
condensates are no n- vanishing. It turns out, however, that the true minimum of the effective 
potential is at a point where one of the condensates vanishes except at the particular point 
5 = (to be discussed further below) where one is at a first order phase transition so that 
phase coexistence can occur. Depending on the sign of 5 we therefore have very different 
behavior at zero temperature and chemical potential- namely either chiral condensates or 
Cooper-pairs being present. 

The phase structure of the limiting theories which have only one coupling constant is 
easy to determine analytically, and we include for completeness a discussion of these two 
particular cases which benchmark our numerical study of the general case. Performing the 
integration in our expression for the renormalized effective potential numerically, we then 
map out the phase diagram for the more general two coupling constant function of 

5. We find a regime of (positive) 5 which remarkably mimics the phase structure of two flavor 
QCD described above. It has a tricritical point as well as a first order phase transition as a 
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function of chemical potential from the chirally broken phase into the superconducting phase. 
We determine the tricritical point and the critical temperature at which the superconducting 
phase transition occurs as a function of 6. As we decrease the magnitude of 6, first the 
regime of first order phase transition from the chiral phase disappears, and then at the 
special point 6 = the possibility of a chiral symmetry broken phase totally disappears so 
that when 6 < one only has the possibility for a superconducting broken symmetry mode 
at low temperatures. 



II. GENERAL CONSIDERATIONS 



We consider the most general Lagrangian with quartic fermion couplings, possessing 
0{N) flavor symmetry and discrete chiral symmetry. 

+ 2G2(^W^g^(j))(^W^^^{i)) _ ^^t«^«. (2.1) 

The flavor indices, summed on from 1 to A^, have been explicitly indicated. The flrst quartic 
term is the usual Gross- Neveu interaction, whereas the second such term, which differs in the 
arrangement of its flavor indices, induces the pairing force to leading order in This term 
is possible because we demand only 0{N) symmetry as opposed to the SU{N) symmetry 
of the original Gross-Neveu model. In the flnal term, /i is the chemical potential. 

Strictly speaking, a (ipip) condensate cannot form, because it breaks the U{1) of Fermion 
number and hence violates Coleman's theorem [|13|. Similarly, (ipip) as well as (ipip) con- 
densates cannot exist at flnite temperature in one spatial dimension because of the Mermin- 
Wagner theorem |TB|. Nevertheless, it is meaningful to study the formation of such conden- 



sates to leading order in as explained in ref. ^4 



Our conventions are: 7° = cri; 7^ = —"1(72', 75 = c"3- The pairing term, proportional to 
G^, may then be rewritten: 

2G2^»75^(^'^^«75^(^'^ = -G2[6,/3^«^V?'^]K5^i^'V?] . (2.2) 



Following standard techniques [8] we add the following terms involving auxiliary flelds m, 
B\ and B: 

A £ = -^[m + g^i^f - i^(5t _ Gh^pi,^^Hf){B + G\,i,^^^ i^^f) . (2.3) 



4 



This addition to L will not affect the dynamics. In = £ + A£, the terms quartic in 
fermion fields cancel, and we have 



a = ^{t t-m- /i7°)V^ -—.-^ + " Bh^.^j^^^^ . (2.4) 



i2 

2^2 ^2 



We integrate out ip and t/;^ to obtain the effective action depending on the auxiliary fields 
m, B and i?^: 

r,^;(m,S,5t) = I d'xi-^ - ^) - '-Tr In A -'-Trlnil + M^A')-' a, a,] (2.5) 

where we have subtracted a constant (independent of the auxiliary fields) and have defined: 
^ = 7°(^ i/~ ^ ~ /^7°) = "^^0 + ^cTs^^z — f^ — "^cTi (2.6) 

so that = —ido — ia^dx — ^ — mai. 

Since we are looking for a vacuum solution, we have assumed in ( |2.5| ) that B, B'^ and m 



are constants and have set = AB'^ B. The trace on flavor indices will give a factor A^. 
The large-A^ limit is achieved by setting g'^N = A, and G^N = k/A, and letting iV — >■ oo 
with A and k fixed. We define the effective potential Veff via 

Te// = -N{ J d'x)V,ff (2.7) 



and we therefore have 

V^ffim, M) = ^ + + KV';(m, M) , (2.8) 



m2 M2 ,^(1) 



with V}j}{m,M) = |[tr In(A^A)^^ + trln(l + M2(A^)- Va^- Vs)^^], where now the trace 
is only over the spinor indices. 

We next generate the local extrema of Veff by solving 

dVeff _ dVeff 



dm? dM'^ 
.(1) 



. (2.9) 



We evaluate the matrix products in V^j^j in momentum space, with — > ik^. The traces 
can be done with the help of 

tr[ K ] = (2.10) 

Vo + V-a V?-V^ 
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for any Vq, V. After some manipulation, equations ( |2.9| ) become 



1 _ 9\/>; _. r (Pk \kl -kl + fj.^ + M' 



m 



2 A dm^ J (27r)2 D 
1 _ dV}jj _. r (Pk [kl -kj-fi'~M^ + m' 



K 9M2 J (27r)2 D 



(2.11) 



(2.12) 



where 



D=[k^- kf -M' + m'- /iT - 4KA;^ + /i'A;f - m'kf] (2.13) 



In this expression, ko is shorthand for ko+iesgnko, where e ^ 0+. This prescription correctly 
implements the role of fi as the chemical potential. 

The equations can be reduced further by doing the ko integral. Let us define k± = 
y/bi ± 262, where bi = M"^ + m'^ + fi"^ + kf , and 62 = [M'^m'^ + fi'^{kl + m'^)]^ . Then evaluating 
the fco integral by contour methods, taking proper account of the ie prescription mentioned 
above, we find 

1 _ 1 /-A 1 1 {M^ + fi"^) 1 1 

2A Stt i-A ^ k+ k_ ^M^m^ + i2'^{kf + m'^) K ^- 

and 

rffci— + — + , • 2.15 



K Stt J-a A;_i_ k. 



The kl integrals are logarithmically divergent and we have regularized them by imposing a 
cutoff A. This will be absorbed in the renormalization process to be described in the next 
section. Note, however, that the combination ^ = ^ — jx given by a convergent integral. 
This fact will ultimately lead to the renormalization-scale independent constant mentioned 
in the introduction. 

We observe from the form of equations ( [2.11 ) and (2.12) that the function Vg^j- can be 



reconstructed by integrating with respect to and in the expressions for 2^ and K 
This will determine Vgjj. up to a single constant v}jj-{0,0), which can be chosen arbitrarily 
without affecting any physical quantity. Explicitly performing this integration we obtain for 
the unrenormalized determinant correction to the effective potential 

V^^\m,M) = dki\k+ + k] (2.16) 
2tt Jo 
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To generalize this discussion to the case of non-zero temperature, one returns to eqns. 
( p.ll| ) and ( p.l2| ), and one continues to Euchdean space via the replacement —ik^ with 



A;4 now considered real. The statistical-mechanical partition function is obtained from the 
Euclidean zero temperature path integral by integrating over a finite regime in imaginary 
time T = it from to /5 = Because of the cyclic property of the trace, the Fermion 
Green's functions are anti-periodic in r and one has the replacement 

dh^'^E (2.17) 

P n 

where the antiperiodicity gives the Matsubara frequencies: 

(2n + l)7r 

uJn = k,„= '- ^ > (2.18) 

To do the sum over the Matsubara frequencies, one uses the calculus of residues to obtain 
the identity: 

|5:/(^a;„) = -^tanh^ Resf{z,) (2.19) 

where Zg are the poles of ji^z) in z in the complex plane; Resf{zs) is the residue of f{z) at 
Zg and we have assumed the function f{z) falls off at least as fast as 1/|2;|-'^"'"'^ for large \ z \. 
It will be convenient to use: 

tanh^^^ = 1 — 2nf{zs) 

where 

1 

nf{z) 



e^^ + 1 

is the usual Fermi-Dirac distribution function. 

Rotating equations (2.11) and (2.12) into Euclidean space as described above, we get: 

1 _ dV^ff _ /■ <Pk l-kj -kj + fi^ + M^-m' 



2A J (27r)2 D 

1 dV}jj f (Pk [-kl - kl - fi^ - + 



K dM^ J (27r)2 D 
where now the integral on k^ is defined by eq. (|2.17| ). 



where D = [—k\ — k\ — + m? — /i^]^ — A[—m'^kl + fikf — m?kf\. There is no longer any 
need for an ie in the definition of /C4. Performing the sums over the Matsubara frequencies 
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we obtain the unrenormalized form of the equations which are given by the same expression 
as the zero temperature ones found earher, with the replacements: 

-L^^(l-2n,(M) 

f -^(l-2nKM) (2.20) 

As before we can integrate this to get the determinant correction to the effective potential 
which in unrenormalized form is: 

V^^\m,M) ^ rciA;i[A;+ + A;_ + ^ln(l + e-^^+) + |ln(l + e-^'=-)] (2.21) 
2tt Jo p p 



III. THE CASE /X = T = 



Renormalization of the effective potential is best discussed in the context of the zero 
temperature and density sector of the theory where we can define the renormalized coupling 
constant in terms of the physical scattering of Fermions at a particular momentum scale. 
This vacuum sector is interesting in its own right and we shall be able, by analytic means, 
to derive the result that depending on a parameter 5 related to the relative strengths of the 
two couplings the theory will be in one or another broken phase and only in a mixed phase 
when 5 = 0. Setting /i = T = we obtain 

dV^ff 1 /-A 1 ,^ M. 



dm? 



-— / dh[{l + —) I + 1 , 3.1 

47rio m + + Mf ^ ^kl + (m - Mf] 

_ 1 dh[{l + ^) , ^ + (1 - ^) , ^ (3.2) 



47rio M ^fc? + (m + M)2 M ^fcf + (m - Mf 

which is solved by 



V^^\m,M)^-^ 1^ dhlJkj + (M + m)2 + Jkf + (M - - 2h] . (3.3) 
27r Jo 

This can be integrated to give the unrenormalized effective potential: 

1 ?A 2A 

- — [(M + mf ln(— ) + (M - mf ln(— -)] . (3.4) 
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We renormalize by demanding that the renormahzed couphngs kr and \r satisfy 

OBOB^ m=mo Kr 



and 



I - ^ (3 6) 

— ^ |m=Mo— t— • 1,0. Dj 

am^ m=mo Ar 

Here M = Mq, m = mo designates an arbitrary renormahzation point on which the coupUngs 
will depend. Using these conditions to solve for A and k in terms of Xr and kr yields the 
renormahzed form of the effective potential: 

Veff ^m^[a + — In \ \] + M^[b + —In 



Att 7o 47r 7o 

1 , M + m 

— mM in — 

27r M — m 



+ —mM In I — I (3.7) 



where a and b are the following constants: 

_ J 3_ 

2A^~4^ 

^' = --- + -^1-1^^^! (3-8) 
27r SttMo ' Mo + mo 

and 7o =| Mq — mg |. 

Note that the renormahzation we have just performed at fj, — T — is also sufficient to 
remove all divergences from the effective potential in the more general case of non- vanishing 
chemical potential and temperature. The addition of /i and T will only result in finite 
corrections to the gap equations and therefore to the vacuum values of m and M. We shall 
return to this point in section V. 

For future reference we also want to consider the special renormahzation point relevant for 
the sector where there is chiral symmetry breaking but no Cooper-pair gap when /x = T = 0. 
That is we will choose our renormahzation point to be the minimum of the potential which 
occurs in that case at 

m = mF; M = 0. (3.9) 
For the choice mo — mp, Mq — , a remains the same but b reduces to 

6= — -A (3.10) 
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The renormalized coupling Ai? takes on the particular value vr and Vejf simplifies to 

Veff = M^i + ^ n . -1 + —mM n — 3.11 

" \r Air' Air ^ ' m| ' ' 2n ^ M - m ^ ^ ' 



Here we want to point out that the quantity 

c- , 1111 , , 

5 = b-a = — = 3.12 

so that S is the same number before and after renormalization. 

The gap equations are properly derived by differentiating Veff with respect to B and m 
and then setting these derivatives to zero. Because Veff depends only on B'^B and m^, it 
will always be possible to have solutions with one of m or i? or perhaps both set to zero. 
Differentiating eq.( |3.7| ) we obtain the gap equations: 

1 1 I -m^ L M , , M-m , ^ 

m\2a H \ In -] In = (3.13) 

and 

, .r, 1 1,1 M'^ — m^L m , , M — m , 

M6+ — + — In^ ^ In =0. 3.14 

^ Att Att 7o ^ 47r ' M + m ' ^ ^ 

The solutions m = m* and M = M* will give us the local extrema of Veff. The first of 
these equations is an identity if m = 0, and the second if M = 0. Also the values m* and 
M* that solve these equations are physical parameters that must be independent of the 
renormalization scale 79. Thus these equations tell us how a and b individually run with 70. 
We note, however, that if we solve for the combination 

S = b-a = —[- ^] In I — I (3.15) 

Att^ m*M* ' ' M* + m* 

the scale 70 drops out. Therefore 5 is a true physical parameter in the theory; we shall see in 
the next section that its value controls which of the two condensates m and M can exist. In 
the particular case where the minimum of the potential occurs when m* = mp and M* = 
we have the simple result: 

= (3.16) 

2Xr Kr 271 
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IV. ANALYSIS OF THE GAP EQUATIONS 



It will be useful in the following to note that, at a solution of the gap equations (3.12) 
and (3.13), the effective potential takes the simple form 

Veffim, M) = -^{m' + M') . (4.1) 

Our goal is to analyze all the solutions of the gap equations and to find the one that produces 
the global minimum of Vg//. This will then represent the true vacuum of the theory. 

There are four types of solution to ( |3.13D and ( p. 141 ). The first is simply to set m = M = 



0, leading of course to = 0. Clearly, from ([4. 1|) we see that if any other solution exists, 
V = cannot be the minimum of V. The second and third types are obtained by setting 
M = 0, m 7^ and m = 0, M ^ respectively. If M = 0, then from ( |3.13|) , we have 



= 70 e-(^+^™) (4.2) 



so 



Vo{m, M = 0) = _^e-(i+^™) (4.3) 
An 



(we shall use Vq to denote values of Vg// at solutions of the gap eqn.). Likewise, if m 
0, M ^ 0, then from (|37[^ ) 



= 70 e-(^+^"') (4.4) 



Vo(m = 0, M) = e-(^+^"') . (4.5) 

4:71 

Thus we see that 

Vo{m = 0, M) < Vo{m, M = 0) if 5 < (4.6) 

and 

Vo{m, M = 0) < Vo{m = 0, M) if 5 > . (4.7) 

The fourth case is when both m and M are non-vanishing. The analysis of this case is 
presented in the Appendix where it is shown that the solution with non-vanishing m and 
M always has V^ff intermediate between the values of V^ff associated with the two cases 
where one or the other of the condensates vanish. 
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We conclude that the global minimum of Vejf has M = 0, m 7^ if 5 > 0, and m = 
0,M ^ if (5 < 0. 

As we shall find later, the special point 5 = is the limit point of the line in /i, T space 
where there is a first order phase transition from the phase with chiral symmetry breaking 
to the phase where there is only superconductivity. 

V. RENORMALIZED EFFECTIVE POTENTIAL 



From eqn. (|2.21|) we can see that the corrections due to non-vanishing temperature and 
density do not affect the ultraviolet behavior of the integrand in the ki integral defining 
V^-^\ Therefore, the renormalization that we have performed at /i = T = in section III 
suffices to remove the ultraviolet divergences from the effective potential, and will allow us 
to send the cutoff to infinity. It is perhaps worth recording the complete result explicitly. 
We find, from eqns. (|3.5| ) and ( p.6| ), that 

^ = « + i + ^ (5.1) 
i = (,+ J_ + X (5.2) 
where a and b are defined by eqn. (^.81), and X is a divergent integral given by 



1 1 1 

X = — dU , + 



Air Jo 



7A;2 + (mo + Mo)2 y^fcf + (mo - Mq)^ 



1 2A 

— [In ( )] + terms which vanish as A ^ 00 . (5.3) 

271" a/TO 



Thus the full renormalized effective potential may be written 



1 f°° 2 

Veff = aim^ + aaM^ / dfcJfc, + A;„ + - In (1 + e"^^'^ 

27r Jo P 



+ |ln (1 + 6-'^'=-) 

~2ki-( + ^ = . 5.4 

2 ^kl + (mo + Mo)2 ^kf + (mo - M,)^ 



where Oi = ^(1 + Ana) and 0^2 = ^(1 + Anb). If ai < a2, then at /i = T = the vacuum 
has m^ = = 7oe~^'^"^ and = 0. Here m^T^ is the dynamically generated Fermion 
mass. It is convenient to choose the renormalization scale so that m|. = 70. This entails 
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setting «! = 0. Furthermore, we are free to choose Mq = 0, so that mo = mp. Then Vejf 
takes the form 

Veff = asM^ - — / dh[k+ + k_ + -(In (1 + 6'^''+) + In (1 + e"^"'-)) - 2h 

ITI Jo fj 

{m' + M') / ] ■ (5.5) 

'k( + nip 



Here 

a2 = S = — 77- > 0; Xr = 11, 

as described above in Sec. III. 

It is this branch of the theory that we are interested in as a model for QCD, since QCD 
at zero temperature has a chiral condensate, but does not have a Cooper-pair gap. We 



observe that if we set M = in this expression, we obtain, with E = \^W+m?, 
K//(m^T,/z) = — [In 

2 dk 



[In (1 + e-'3(^+^)) + ln(l + e'^^^-'^))] (5.6) 

P JO 271 



which is the effective potential for the Gross-Neveu model in agreement with refs [0 [jT6 



We will use the analytic information already known about the G-N model as a benchmark 
for our numerical work below. 

In the opposite case a2 < ai, we have, in the yU = T = vacuum, = and = 
= 7oe~^'^"2, -where A is the dynamically generated gap. So we choose 02 = 0, ai > 0, 
and rriQ = 0, A^ = 70 = Mq. The effective potential becomes 

K// = am^ - — / dk^[k+ + k_ + -(In (1 + e-^'=+) + In (1 + e"^^-)) - 2k^ 
27r JO p 

-{m' + M') ^ / J . (5.7) 
kl + A2 



For this case, by choosing 0:2 = we obtain: 

kr = Att, -6 = ai = 77- > 0- 

ZAji Zn 

When = 0, this expression gives us the effective potential at finite temperature for 
the pure Cooper-pairing model considered in |^. Explicitly we have 

Note that it is independent of the chemical potential, as was the case at T = 0. 
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VI. PHASE STRUCTURE OF THE CLASS OF MODELS 



A. Cooper Pair Model 

The pure Cooper pair model has the property that the chemical potential is irrelevant 
and can be transformed away. The form of the effective potential is exactly the same as that 
for the Gross Neveu model at zero chemical potential with M replacing m and the gap A 
replacing mp. Thus in leading order large-N there is a second order phase transition to the 
unbroken mode at a critical temperature which can be determined by the high temperature 
expansion. For T » M we can expand the integral in eq. (5.8) to obtain 



K// = — [In (^)-7], (6.1) 

where 7 is Euler's constant. The minimum of this function occurs at M = 0, which means 
that the condensate vanishes for large T, as expected. The critical temperature is that 
temperature for which 

In (^) - 7 = ^ T, = -e\ (6.2) 
The same critical temperature was obtained in another variant of the Gross- Neveu model 



which had a superconducting phase [0, so that this temperature seems ubiquitous in 4- 
Fermi models in 1+1 dimensions. 



B. Gross-Neveu sector 

As is well known, the Gross Neveu Model has spontaneous symmetry breaking at zero 
chemical potential and temperature. At zero temperature, the symmetry is restored at finite 
chemical potential at a critical value of yU [ITD] |TT|. This transition is first order. At zero 



chemical potential the system undergoes a second order phase transition to the unbroken 
symmetry phase as we increase the temperature. Thus at some point in the phase diagram 
there is a tricritical point. For this model we have performed both high and low temperature 
expansions of the leading order in 1/A^ potential which is given by: 

2 2 
2 „ V m „ m 



Veffim ,T,fi) = —[In ^-1] 

-I [In (1 + e-'^^^^^)) + ln(l + e'^^^"'^))] (6.3) 

p Jo Ztt 
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In the high temperature regime, using methods similar to those used for Bose condensa- 



tion ITSll we obtain: 



^e=^e^[l-^^] (6.6) 



K//K,r,/i) = — (In ^ + 2^(/^ (6.4) 
which leads to the relationship: 

Tc = —-Ml-^-^] (6.5) 

which at /i= gives the same critical temperature as for the Cooper pair model, however 
with mp replacing A. At small one has approximately 

_ 7/i^C(3) . 

In the low temperature regime we want an analytic expression for the effective potential 
which would enable us to determine values of /i and T at which the first order transition 
occurs. At zero T the modification to the effective potential due to the chemical potential is 
only in the region m < fi. The standard low temperature expansion used in Bose Condensa- 
tion [0 unfortunately only gives the finite temperature corrections when fi < m and thus 
is not very relevant to the question we want to answer. To obtain an approximate analytic 
expression valid in the opposite regime m < fi pertinent to the first order phase transition, 
we resort to a crude approximation which captures the relevant physics. That is we make 
an approximation to the Fermi-Dirac distribution function that allows us to perform all the 
integrals. First we rewrite the derivative of the potential in the form: 



dV m ^ m? m f^dk, ^ E + ji ^ E — ii, , ^ 

— = — ln^ + — / — 2-tanh — - tanh — , 6.7 
dm 2tt ml tt Jo E^ 2T 2T ^' ^ ^ 

where E = y/k'^ + m^. and then replace the function tanh(^^^) using the straight line 
interpolation: 

tanh(a;) ^ {1 if a; > 2; - 1 if x < 2; a; if |a;| < 2}. (6.8) 

This has the correct behavior as T ^ and captures the physics of the broadening of the 
Fermi surface. At T = 0, the effect of the chemical potential is the most dramatic. Since 
in that limit tanhx = e(x), we get immediately that 

dV m m^ m fVi^^-^dk^, 

^^Zl2+— L — e(/i-m) 



dm 2ti mp n Jo E 

,2 ^ , 1 ^2 



m , m m ^ , ^ r, , I ^ , 1 , m ^ 

— In — + -0(/x - m){ln(l + m^ j fj?) - - In — }, (6.9) 

/TT mp vr * Z fi^ 
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This can be integrated to give the result that for m < jj, the effective potential is given by: 



K// = ^{m\2 ln[ ^ + ^^' - 1) - 2^iJii^ - m^} + (6.10) 

47r nip 

whereas, for m > the effective potential is equal to its = value, namely 

2 2 
Tfi Tfl 

^e//= 4^[ln ^-1]+C(/.) (6.11) 

The arbitrary integration constant can be eliminated by choosing Veff{m = 0) = 0, which 
yields 

C(/i) = ^ (6.12) 



For all T we can use the approximation in Eq. ( |6.8| ) to perform all the integrals explicitly. 
Doing this, we obtain an approximation to the exact phase structure in the regime where 
there is a first order phase transition as is shown in Fig. 1. In that figure we also include the 
high temperature analytic result. Our analytic calculation gives us an approximate value 
for the tricritical point which separates the regime between the first and second order phase 
transitions: — = .661 — = .31 compared to the "exact" numerical result as for example 
found in Ref. 

= .608, — = .318. (6.13) 

mp nip 



C. Full Phase Structure 

The phase structure is quite different depending on whether we choose the case 6 > 
which has chiral symmetry breaking in the vacuum, or 5 < where there is Cooper pair 
formation in the vacuum. In the regime where S > 0, the phases of this model are quite 
similar to QCD as shown in figures 2 with the value of a2 = S = In the vacuum there 
is chiral symmetry breakdown. As we increase the chemical potential at low temperatures 
there is a first order phase transition into a phase with Cooper pairs. At and near the phase 
transition line there can be coexistence of the two separate phases, one with Cooper pairs 
and one with a chiral condensate which breaks chiral symmetry. For the range 

. 1-13097 , 

6 > — = 6, 6.14 

the theory will have a tricritical point at the value given by Eq. ( |6.13| ), so that the regime 
where there is chiral symmetry breakdown will, for chemical potentials below the tricritical 
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value undergo a second order phase transition at large temperatures. For values of the 
chemical potential between the tricritical value and the value for the first order transition 
to the superconducting phase (determined below), the phase transition from the chirally 
broken mode to the unbroken mode will be first order at large temperatures. 

As we move to higher values of the chemical potential, for sufficiently low temperatures 
the system exists in a superconducting phase with nonzero gap. As we increase the tem- 
perature at fixed large chemical potential, the system undergoes a second order transition 
into the unbroken mode, with the critical temperature depending only on 6 and not fi. This 
dependence is displayed in figure 3. Tc reaches the tricritical value Tc/mp = .318 when 
S = Sc- Figure 2 is in the regime where 6 > 6c so that it displays a tricritical point. For 
values of S < 6c, the chirally broken phase only can be restored via a second order phase 
transition. This case is illustrated in fig. 4 which is for 5 = ^ < 6c- In between the chirally 
broken and superconducting phases is a coexistence curve. The intersection of this curve 
with the line T = can be determined as a function of 6 which we will shall do below. 
The existence of two phases having the same energy is shown in the 3D plot of the effective 
potential as a function of m, M in fig. 5 and in the two dimensional slices of this figure shown 
in figs. 6 and 7. The particular case displayed is for a2 = 6 = All the plots are for 
T = 0.02, /i = 0.56 which is numerically determined to lie along the first order line sepa- 
rating the chiral condensation phase from the Cooper condensation phase. As 6 approaches 
zero, the coexistence curve approaches /i = and after that one no longer has a phase with 
chiral symmetry breakdown. 

Infinitesimally to either side of the coexistence curve we are at two separate minima of 
the potential. In each phase the minimum takes place with the value of the other condensate 
mass equal to zero. Thus the condition defining the coexistence curve is 

V[m = mF,M = 0]=V[m = 0,M = M*] (6.15) 

The value of M* is chosen to minimize V^[0, M] for a given value of /i, 6. Recall that at T = 
the effective potential is given by: 

Veff = 6 - — r dk[\ + m2 + M2 + ^2 + 2 Jm?M'^ + (F + m2) /^2 
27r JO * 

+ \jk^ + m2 + M2 + /i2 _ 2 Y/m2 + (P + m^) fi^ - 2k - ] + C(/i) (6.16) 

~\~ k"^ 

We notice that when M = 0, this potential becomes that of the Gross-Neveu model. 
Thus if we choose Ve//(m = 0, M = 0) = 0, then again C(yu) = as in eq. ( |6.12|) . 
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At T = it is possible to analytically determine the value of the chemical potential as a 
function of S as well the value of M at the minimum. On the left hand side of the coexistence 
we need to evaluate the GN effective potential in the regime where m — mp > since the 
phase transition to the superconducting phase always occurs in that regime. Thus we have 

K//(m = m^,0) = -^ + ^ (6.17) 

On the right hand side we need to evaluate the zero temperature effective potential for 
m = 0, M = M* . We have on the Cooper condensation side, 

Veffim = Q,M) = 5M^ + (In -1) (6.18) 

47r rrip 

The quantity M* is determined by that value of M that minimizes this function, namely 

dM 47r ^ ' 

or 

= mle-^''^ (6.20) 

Inserting this value into the equation equating the value of the potential on both sides of 
the phase transition we then obtain the critical value of the chemical potential 

2 1 _ -47r5 



VII. CONCLUSIONS 

In this paper we have analyzed a 1+1 dimensional model possessing 0{N) flavor sym- 
metry and discrete chiral symmetry, and have found a phase structure remarkably similar to 
that conjectured for 2-flavor QCD. We have derived the general forms for the effective poten- 
tial in leading order in j^. We have analyzed the case /j, = T — analytically, showing how 
the phase structure is governed by the renormalization group invariant S. For /j, = T = this 
structure is remarkably symmetric in the two condensates m and M. We have performed 
careful numerical analysis of the integrals involved in the determination of the effective po- 
tential and have determined the dependence of Ve// on the parameters 6, fi and T. What we 
have found is that when there is chiral symmetry breakdown in the vacuum sector { 6 > 0), 
there are at least three different regions. In the low temperature regime, as we increase the 
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chemical potential there is a first order phase transition to a regime which has a Cooper 
pair gap (superconductivity) but no chiral symmetry breakdown. Along and near the phase 
transition line, there is a regime where the two phases coexist like ice and water. At T = 
we explicitly determine the value of the chemical potential at which this occurs and also 
the value of the Cooper pair gap as a function of 5. At high enough temperatures both 
symmetries are restored. In particular US > Sc— then there is also a tricritical point 

so that depending on the value of /i the phase transition out of the chirally broken phase will 
be either first or second order. We illustrated the phase structure of this model by showing 
the phase diagram of this model as a function of temperature and chemical potential for 
representative values of S. We also plotted the effective potential at a representative place 
where there is phase coexistence. In the opposite case S < one finds that in the vacuum 
sector the theory has a Cooper pair gap but no chiral symmetry breaking. In that case the 
theory has a transition at high temperatures to the unbroken mode where the gap goes to 
zero. 

Using this toy model we intend to study how the phase transition from the high temper- 
ature to low temperature regime proceeds in time during an expansion of an initial Lorentz 
contracted disc of quark matter starting from various initial conditions related to differ- 
ent points on this phase diagram. Wc hope to determine how various correlation functions 
depend on the initial conditions of a scattering experiment, assuming that it produces an 
initial state in local chemical and thermal equilibrium somewhere on the phase diagram we 
obtained in this paper. 
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APPENDIX A 



In this Appendix we give the details for determining that the relative minimum which 
has both condensates is always between the two minima which have only one condensate. 
When both m and M are non-vanishing, it is then convenient to define p = — and to 
combine the gap equations in the form 



and 



1 (l-p2) 



Att p P + 1 



In I I . (7.2) 



Both these equations are even in p, so we may take p > for convenience. Eqn. ( [7.2|) tells 
us immediately that if5<0,0<p<l, and if 5 > 0, p > 1. Furthermore, the r.h.s. of 
( [7. 21 ) is bounded between and Hence we conclude: If | 5 |> ^ there is no solution 
with both m and M non- vanishing. If | 5 |< there is such a solution, with the property 
that m > M if (5 > and M < m if (5 < 0. 

It remains to decide whether Vo{m,M) can be the global minimum. To this end, it is 
convenient to re-express the gap equations once more in the following form: 

- (1 + Ana) = In ' ~ ^ ' - ln{ | ^ (7.3) 

7o P + 1 



and 



|p^ — 1| ,r,P— 1,1^ 
4nb) = In ^— ^ - ln{ | |p} . (7.4) 

7o P + 1 



From these, making use of eqs. ( [4.1|) , ([4.3|) and ( [4.5|) , we immediately obtain 



and 



where 



Voim = 0, M) = e-(i+^-'') = (7i(p)Vo(m, M) (7.5) 

4:71 



Voim, M = 0) = e-(i+4™) ^ ^2(p)K)(m, M) (7.6) 



(l + p)^+) |l-p|^-^ 

^i(^) = TT? ^ ^ 
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and 



(7.8) 



Eq. (7.5) is the relevant comparison li < 5 < Q and < p < 1, whereas eq. (7.6) is 
relevant for < 5 < ^ and p> 1. 

We observe, however, that g2{^) = 9i{p): so both cases reduce to the following: if we 
can show that gi{p) > 1 in the range < p < 1, then Vo{m, M) is never the global minimum 
(recall that the VqS are all < 0). On the other hand, if gi{p) < 1 in this range, it will be 
possible to have Vo{m, M) be the global minimum. 

To settle this question, write gi = e^, with 



In the range of interest, p^ < p, so the r.h.s. is a sum of positive terms. Hence h{p) > and 
9i{p) > 1- 

We conclude that the global minimum of Vg// has iV/ = 0, m 7^ if 5 > 0, and m = 
0, M 7^ if 5 < 0. 



Mp) = (1 + -) In (l + p) + (l--) In (1-p)- In (1 + p') 



P P 



(7.9) 
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FIG. 1. Phase diagram for the Gross- Neveu model. Partial hnes are results of 
the High and (approximate) Low Temperature expansions. Continuous line is 
the numerical result. 
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FIG. 2. Phase Structure at S — Tricritical point is at T/mp — .318, fx/niF — .608. 
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FIG. 3. Critical Temperature for the Superconducting Phase as a function of 0:2 = 5 
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FIG. 4. Phase Structure a.t 6 = j- 




0.25 



0.75 



FIG. 5. Phase Coexistence Effective Potential at T = 0.02, yU = 0.56,5 = j- 
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FIG. 6. Phase Coexistence Effective Potential as a function of M for m = 
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FIG. 7. Phase Coexistence Effective Potential as a function of m for M = 
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